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A quantum processor is a programmable quantum circuit in which both the data and the program, 
which specifies the operation that is carried out on the data, are quantum states. We study the 
situation in which we want to use such a processor to approximate a set of unitary operators to a 
specified level of precision. We measure how well an operation is performed by the process fidelity 
between the desired operation and the operation produced by the processor. We show how to find 
the program for a given processor that produces the best approximation of a particular unitary 
operation. We also place bounds on the dimension of the program space that is necessary to 
approximate a set of unitary operators to a specified level of precision. 
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I. INTRODUCTION 

Quantum circuits are typically designed to perform one 
function, for example, teleportation or cloning. It is use- 
ful to have circuits that are more flexible and can perform 
a variety of functions. The operation that the circuit per- 
forms can be determined either by classically setting the 
values of some parameters, for example the rotation an- 
gle in a one-qubit rotation gate, or it can be determined 
quantum mechanically, where a quantum system serves 
as a program to tell the circuit what to do. The second 
method has the advantage that the program could be 
a result of a previous stage of a quantum computation, 
which would allow one stage of a computation to control 
a subsequent one. 

A programmable quantum circuit (quantum processor) 
has two inputs, the data register and the program regis- 
ter. The data register is in the state on which we want 
to perform an operation, and the program state specifies 
the operation. An example is the C-NOT gate in which 
the control qubit is the program and the target qubit is 
the data. If the control qubit is in the state |0), nothing 
is done to the data state, and if it is in the state 
the operation Ux (bit flip) is applied to the data state. 
The superposition state q;|0) -f /3|1) causes the completely 
positive quantum map 

T{p)^\a\^p+\l3\^cTxP<Jx, (1.1) 

to be applied to the data state, p = \%l]){'4)\. 

Suppose that we have a set of N unitary operators that 
we want to be able to implement on the data qubit with 
a programmable quantum circuit. Nielsen and Chuang 
showed that this requires a program space of at least 
N dimensions [1]. This follows from the fact that the 
program states corresponding to any two of the unitary 
operators must be orthogonal. If one wants to be able to 
realize a large (or infinite) number of unitary operations 
with a program space of fixed dimension, one has two 
possible options: One option is to make the processor 



probabilistic, that is a measurement is performed at the 
program output, and if the correct result is obtained, the 
desired operation has been performed on the data [1]- 
[5] . The probability of obtaining the proper measurement 
outcome will, in general, be less than one so that the 
processor succeeds with only a certain probability. The 
second option is to make the processor an approximate 
one. That is, each of the operations is not performed 
exactly, but only up to some level of approximation. It 
is this type of processors that we wish to discuss here. 

Approximate processors have been discussed by Vlasov 
[2] and by Vidal and Cirac [3] . Vlasov considered a clas- 
sically programmable processor, while Vidal and Cirac 
considered one whose programs are arbitrary quantum 
states. They made some rough estimates of the resources 
required for a processor to be able to program a set 
of unitary operators to a specified level of precision [3]. 
In Refs.[6]-[8] approximate programmable quantum mea- 
surement devices have been studied. These devices real- 
ize certain classes of POVM's up to some level of approxi- 
mation, and which POVM they perform is determined by 
a program state. It was shown by D'Ariano and Perinotti 
that for programmable measurement devices the number 
of dimensions of the program space is a polynomial func- 
tion of the reciprocal of the desired accuracy [8] . 

Physical limitations in real systems lead to an addi- 
tional reason to study approximate processors. That is, 
ideal devices in theory become approximate ones in prac- 
tice. For example, if one wants to perform a rotation on a 
qubit that is a two-level atom, one applies a classical field 
to the atom. Real fields, however, consist of photons, and 
this and energy constraints on the field place limits on 
the accuracy of the rotation that can be achieved [9]-[12]. 
Conservation laws can also place limits on the accuracy 
of quantum operations [13]. 

Other types of processors have been explored. In par- 
ticular, processors that evaluate the expectation value 
of an arbitrary operator have been proposed [14, 15]. In 
these processors, the data is the state in which the expec- 
tation value is to be evaluated, and the program specifies 
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the operator. 

In this paper we shaU discuss processors that approxi- 
mate sets of unitary operators. Wc shall show, for a given 
processor, how to select an optimal program vector to ap- 
proximate a particular unitary operator. In addition we 
shall give a lower bound on the number of dimensions 
the program space must have to approximate a set of 
unitary operators to a given level of accuracy. In the last 
part of the paper wc will address the question of opti- 
mal programmability, i.e. which processor is the best in 
approximating all channels. 



II. OPTIMAL PROGRAM STATES 

We now consider a processor that acts on the Hilbert 
space Ti. = Tid (S) Hp, where Hd is the data Hilbert space 
and Hp is the program Hilbert space. Let us denote 
the dimension of Hd by D and that of Hp by A''. The 
processor itself is represented by a unitary operator G, 
which acts on H. The action of the processor on the 
input state \ip)d\'^)p is given by [16] 

N 

G(|V)d®|S))=5^A,-(S)|VMj)p, (2.1) 
i=i 

where {|j)p|j = 1, . . . N} is an orthonormal basis of Hp. 
The operators Aj (S) are expressed in terms of the oper- 
ators Ajk, where G is expressed as 



TV 

G= ^jk^\j}pp{k\. 
These operators obey the relations 

N 

j,k=i 

N 

k=l 



(2.2) 



(2.3) 



where Id is the identity operator on Hd- The operator 
Aj{E) is given by 



N 



k=l 



from which it follows that 

AT 



^A](S)A,-(S)=/rf. 



(2.4) 



(2.5) 



Wc now need to discuss how to measure how close our 
processor comes to achieving a particular unitary oper- 
ation. We shall use, what has been called by Gilchrist, 
et al., the process fidelity [18], which was originally pro- 
posed by Raginsky [19]. It is defined as follows. Let 



Ti and T2 be two completely positive maps, which map 
operators on the Hilbert space IC onto operators on the 
same space. We shall assume that the dimension of IC is 
finite and equal to D. The Jamiolkowski isomorphism al- 
lows us to associate a density matrix on /C (8) /C with each 
of these maps. Define the maximally entangled state 



(2.6) 



where = 1,...N} is an orthonormal basis of /C. 

For each map Tj, define the density matrix pj to be 



p, = (I®T,)(|$)($|), 



(2.7) 



for j = 1,2, where T is the identity map. The process 
fidelity is defined as 



Fproc{Tl,T2) — 



Tr 



(2.8) 



The process fidelity has a number of useful properties 
that are discussed in Refs. [19] and [18], one of which 
is the fact that it is symmetric, i.e. -Fproc(7i, T2) = 
Fproc{T2, Ti). 

We are going to be interested in the case in which one 

of the maps is unitary. In particular, let us assume that 

Ti(p) = UpU~^ for some unitary operator U. In this 

1/2 

case we have that pi is a pure state so that = Pi- 
This gives us that 



1 1/2 

D ' 
^ (il|C/-'T2(|il)(i2|)C/|i2) 

jl>j2=l 

(2.9) 

If T2 is the result of the action of a processor, we have 
for a density matrix pd, representing a data state, that 



Up) = J2M^)pdM^)\ 
3=1 



(2.10) 



which gives us, finally, that (we denote the map Ti by 
the operator U) 

mT2) = -^^|Tr(C/-%(-))|'. (2.11) 
i=i 

Using the notation for the Hilbert-Schmidt scalar prod- 
uct {A\B) = TrAtS this can be rewritten in the form 

This fidelity can also be expressed in terms of the op- 
erators Ajk - Defining the matrix 

1 ^ 

Mkrk2 = -^J2^(4k.UmU-'Ajk,), (2.12) 
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we have, from Eq. (2.4), that 



N 

F{U,T^)= ^ p{E\h)pMk,k,p{k2\^)p- (2.13) 

/Ci,fe2=l 

Now consider the following problem. Suppose we are 
given a processor and we wish to find the best pro- 
gram to approximate the unitary operator U, where 
by best we mean the program that maximizes the pro- 
cess fidelity. An examination of Eq. (2.13) shows that 
this can be accomplished by finding the eigenvector of 



\ki){k2\ with the largest eigenvalue. 



and choosing the program vector to be this eigenvector. 
The corresponding fidelity will just be the largest eigen- 
value of M. 

This procedure is particularly simple to carry out when 

the processor is, what was called in Ref. [16], a U pro- 
cessor. This is a processor that is a controUed-U gate. 
Each basis vector \k)p in Hp is associated with a unitary 
operator Uk acting on Hd- That is. if the program state 
is \k)p, then the operator Uk is applied to the data state. 
The operators Ajk for this type of processor are particu- 
larly simple, Ajk = SjkUk, which implies that the matrix 
M is given by 



Mk,k2 



1 

D2 



|Tr(U^UkJ| (5, 



felfe2 ■ 



(2.14) 



Because in this case M is diagonal, we simply find the 
diagonal element that is largest. This is the largest eigen- 
value of M and the maximum value of the fidelity. The 
value of k corresponding to this diagonal element tells us 
which of the basis vectors \k)p is the program that will 
achieve this fidelity. This implies that to best approxi- 
mate a unitary operator ?7 by a U processor, we simply 
find which of the unitary operators that the processor can 
perform perfectly has the largest Hilbert-Schmidt inner 
product with U and perform that operation. Note that 
this prescription does not make use of superpositions of 
the basis states in the processor. 



III. AN EXAMPLE 

Before proceeding with the exploration of the gen- 
eral properties of approximate quantum processors, it 
is useful to analyze the following example. We shall 
consider a processor acting on qubits with an N dimen- 
sional program space spanned by the orthonormal basis 
{|A:)p|fc = 0, ...A'' — 1}. Define the shift operators 
and acting on the program space as E+\k) = \k + l) 
and E-\k) = |fc — 1), where the addition and subtraction 
are modulo N. We also define the program states 



1 



N-l 
fe=0 



-ike 



\k). 



(3.1) 



If = 9m = {2TTm)/N then the state \6m) becomes an 
eigcnstatc of i?+ and E- 



E+\6m) = e^'^-lOm) ; 



E_\err^)=e-'^-\0m). (3.2) 



For the qubit, whose Hilbert space is spanned by the two 
orthonormal vectors |0)d and \l)d, define the operators 
ct(+) and cr(-), where (j'^+'>\())d = |l)d, fT(+)|l)rf = 0, and 
cr(~) = (cr(+))'''. We shall consider a specific realization 
of the U processor defined by the operator G acting on 



G = exp 



(cr(+) (S)E.+a 



(-) 



'E. 



(3.3) 



The fact that G is a U processor can be seen when we 
let G to act on the state \ip)d\0m)p- Here we obtain the 
result 



= G(|V)d |^m)p) 



= exp 



(5) « 



Defining 
U{0) 



exp 



e iV(+)+e'V-)) 



(3.4) 

|V')d O |6'm)p- 

, (3.5) 



we see that we can perform U{9) perfectly when 9 = 
9m, for some m. Suppose, however, we are interested in 
using this processor to approximately perform U{9), for 
6 not equal to any of the 6m- We know what the optimal 
strategy is from the previous section, find the operator 
U{9m) which has the greatest overlap (in the sense of the 
Hilbert-Schmidt inner product) with U{9) and perform 
that operation. Here we are going to examine a strategy, 
which is simpler to implement, but not optimal. We shall 
simply use the state \9)p as a program state. We find that 
this gives us a process fidelity of 



N-l 



= -y 



COs''{9m-0y- 



m=0 



.m^[Ni9m-e)/2] 
sm^[{em-0)/2] ' 



(3.6) 



This sum is an oscillatory function of 9 with a period 
2n/N. The minima of this function are achieved for 9 = 
■K/N+2'Kk/N when the process fidelity takes the minimal 
value Fmtn = l-2/N. 

Let us see how this compares to using the optimal pro- 
gram states. The process fidelity between the operators 
U{9i) and U{92) is given by 



F(f/(^l),C/(^2)) = COs2(^l-^2). 



(3.7) 



If we approximate U{9) by U{9m), where m is chosen so 

that U{9) and U{9m) have the largest Hilbert-Schmidt 
inner product, then the fidelity is bounded below by 

^s-^(^)~>-(^)'- 

Note that in this case the error is of order l/N"^, while in 
the previous case it was of order 1 /N , so there is a cost 
to not using the best program states. 

What we then have is a an approximate processor that 

can be made very accurate by choosing N large enough. 
It achieves an accuracy of order 1/N m. approximating 
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U{d) with the simple program state 16*)^, which is not as 
good as the best accuracy, 1/-/V^, but the approximation 
in none the less a good one for N sufficiently large. Thus, 
we see that a U processor, making use of a simple pro- 
gram, can be quite useful in approximating the action of 
a set of operators labeled by a continuous parameter. 



IV. BOUND ON DIMENSION OF PROGRAM 
SPACE 

We would now like to find a bound on the resources 
required to achieve a given accuracy in approximating a 
set of unitary operators by means of a fixed processor. In 
particular, wc want to sec how the dimension of tlic^ pro- 
gram space grows as the accuracy of the approximation 
increases. 

The Schwartz inequality \{A\B)\ < y/{A\A){B\B) im- 
plies that 



|Tr([/U,(S))| < VD[TriAUE)AjiE))]'/' 



(4.1) 



and, therefore, if the action of our processor with the 
program state is given by the map T, we have that 



F{U,T) = -i^^|Tr(C/tA,.(s))p 

1 N 

< -^Tr{A]{E)AjiE)) = l. (4.2) 

In the last equality we used the normalization property 
of Kraus operators (2.5), i.e. J2j = ^■ 

We begin by assuming that the fidelity is 1 and see- 
ing what this implies about the operators Aj(E). If 
F{U,T) = 1, then, we see from above, that Schwartz 
inequality has to be saturated. This means that the oper- 
ators Aj(E) and U arc colincar, i.e. Aj(E) = 0jU, where 
f3j is a complex number. Furthermore, Eq. (2.5) implies 

SjLi = 1- Now suppose that we have two differ- 
ent unitary operators that can be realized perfectly, Ui 
by the program state and U2 by the program state 

IS2). Therefore, Aj{Ei) = /JyC/i and Aj{E2) = /?2jC/2- 
We then have that 

N N 

J2KjI32jU,-'U2 = ^A](Si)A,(S2) 

= /dp(Si|S2)p, (4.3) 

where wc have used Eqs. (2.4) and (2.3). If Ui 7^ 
1/2, then this equation implies that both p(Si|S2)p and 
SjLi are zero. This result is simply a restatement 

of the Nielsen- Chuang theorem: If two unitary operators 
are realized perfectly by a processor, their program vec- 
tors must be orthogonal. 

Now let us suppose that the processor performs the 
operation U with a fidelity greater than or equal to 1 — 



e, i.e. F{U,T) > 1 — e, where T is specified by Kraus 
operators Aj (S) . Let us express these operators as 



A,{E)^^,U + B,iE), 



(4.4) 



where Tt{U^ BjCE,)) = 0. This decomposition is unique. 
The inequality F{U,T) > 1 — e implies the following con- 
dition on coefficients /3j = ^{U\Aj{E)) 

N N 

1 > F{U, T) = ^ |(C/|A,(S))|2 = ^ > 1 _ e . 



(4.5) 

Tracing both sides of the normalization condition 
"^^j Aj{E)^ AjiE) = I we obtain the inequality 

Tr[i?,(S)ti3,-(S)] = E,(B,(S)|B,(S)) < De. 
Next consider the situation in which our processor can 
approximate two unitary operators, Ui and U2, each with 
a fidelity greater than or equal to 1 — e. In particular, if 
Ti is the map produced by the program state \Ex)p and 
T2 is the map produced by the program state |S2)p, then 
both F{Ui,Ti) and F{U2,T2) are greater than or equal 
to 1 — e. We also have that 



A, (Si) = f3,,Ui+B,,{E,) ; 

A,{E2) = P2jU2+B2j{E2) , 



(4.6) 



where Tv{UlBij{Ei)) = Tr(f/|B2j(S2)) = 0. As in the 
case when the unitary operators were performed per- 
fectly, consider the quantity 



N 



I,{E,\E2) = ^A,(Si)U,(S2) 



(4.7) 



N 



= El/^yf^i + Blj{E,)][(32jU2 + B2jiE2)]. 

Let us evaluate the absolute value of traces of both sides 

D\{E,\E2}\ = \ J2{M^i)\M^2))\ (4.8) 
j 

= I Yllf^iAi^U^\U2) + PljmB2j) 
j 

+f}2j{Bij\U2) + {B,,\B2m 

< \{Ui\U2)\J2 ^hP^o I + + De. 

3 

In the last line we used the formulas 

E,l(5y|52,)l < E,V(Sl,|Sl,)(52,|i?2,) 

^E,(Sy|By)E,(S2,|B2,-) (4.9) 



< 

< De , 



and 



iE,/3r,-(c/iiB2,-)i < T.A0ioWmui){B2j\B2i) 

< D^e. 

(4.10) 
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As a result wc obtain the bound on the inner product 

between two program states 

KHl|S2)|<l|([/l|C/2)||^/3l^/32,■|+2^/i+e. (4.11) 



Next we will estimate the first term. The idea is to use 
Eq.(4.7) and apply both sides to a special vector l'^^) 

that maximizes the quantity 1 — |(^|f/|[72|V')P- Let us 
denote this maximum by 77, i.e. 



r? = max 



(4.12) 



This quantity describes the distinguishability of two uni- 
tary transformations, and a short calculation shows that 
1] < \\Ui - C/2|P- After applying both sides of Eq. (4.7) 
to \ijjri) wc find the components of the resulting vectors 
orthogonal to l-^^) by applying the projection operator 
= I — Itpri) i'tpnl to both sides. The left side vanishes 
and we obtain the equality 



= P,^(E,-4.(Si)U,(S2))|V',) 

= EjPtjP2jP,^UlU2\i>r,) + \c0), 



(4.13) 



where 



N 



+Bl.{E,)B2j{E2)) IV,). 



(4.14) 



We now want to find a bound on ||w||. Using the 
facts that the operator norm is bounded by the Hilbert- 
Schmidt norm, we have that 



N 



|^st(Si)S,.(S2)|| < ^{Bj{E,)\Bj{E,)y 



/2 



{B,{E2)\BjiE2)y/' 
< eD, (4.15) 



and 



N 



N 



\J2PhUlB2ji^2)\\ < J2\0M{E2)\BjiE2)y 



/2 



(4.16) 



Applying these inequalities we have that \\uj\\ < eD + 
2^D. In addition, we find that \\P^U\U2'^rt\\ = y^- 
Therefore, we can conclude 



Defining 




F = min I 1 



(4.17) 



(4.18) 



we have, finally, that 



< ^\iUi\U2)\ + 2Ve + e. (4.19) 

Note that in the case that both operations are carried out 
without error, in which case e = 0, this inequality implies 
that the program vectors must be orthogonal, recovering 
the known result. 

Now suppose that we have M unitary operators that 
we want implemented by a processor so that the pro- 
cess fidelity for each of the operators is greater than or 
equal to 1 — e. How many dimensions must Hp have? 
In order to answer this question, we first find the val- 
ues of Yjfe = {F/D)\{Uj\Uk)\ corresponding to each pair 
of operators in our set, and use these values to find the 
largest set of linearly independent vectors in the set of 
program vectors. Linear independence can be deduced 
from the following result: If {vk\k = 1, . . ■ K} are vec- 
tors of length 1, and K^^/cilffe)! < 1)j then the 
vectors {vk\k = 1, . . . K} are linearly independent [3, 16]. 
Suppose that there is a subset of our operators, with M' 
members, whose pairs have small values of Yj^, and let 
the largest value of Yjk for this subset be Y^ax- Then we 
have for all of the program vectors corresponding to this 
set, that 



<Ymax + 2V~e + e = q{Ymax,e). (4.20) 



Let Kg be the largest integer such that Kg < (l/q) + 
1. What the result we just quoted implies, is that any 
set of vectors whose size is Kg or less, will be linearly 
independent. Therefore, if M' < Kg, then all of the 
program vectors will be linearly independent, and the 
dimension of Tip must be at least M'. If M' > Kg, then 
the dimension of Hp must be at least Kg. This, then, 
is the restriction our result imposes on the dimension of 
the program space. 

As an example, suppose we want to implement the 
operators /, ai,a2, and 0-3 on qubits, where the operators 
(7j, for j = 1, 2, 3, correspond to the usual Pauli matrices. 
For all pairs of these operators wc find that Yjk = 0, and 



q{0, e) = 2\/2i + e. 



(4.21) 



Our bounds then give us that for e < 0.02 the program 
space must have four dimensions, for e < 0.05 it must 
have at least three dimensions, and for e < 0.17 it must 
have at least two dimensions. 



V. ONE-PARAMETER GROUP: TWO 
APPROACHES 

Programmable processors can be exploited to imple- 
ment quantum maps probabilistically. In this case a spe- 
cific measurement on the program state is performed and 
if an a priori defined result is obtained then we know that 
a desired operation has been performed on the data. In 
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other words the specific measurement that is accompa- 
nied by a post-selection induces the desired transforma- 
tion of the data register. As was discussed in [3] a proba- 
bilistic processor without measurement can be used as an 
approximate processor. In this case the transformation 
can be expressed as 



^d^] = PsuccessT[Q] + Perror-J^iol 



(5.1) 



where T is the channel we want to approximate, and 
Psuccess and perror are independent of the input data 
state, p. Due to the concavity of the square root of the 
process fidelity we find that Psuccess < F{£iiT)i i-S- the 
accuracy of the approximation is bounded from below by 
the probability of success. 

Here we want to compare the performance of a prob- 
abilistic processor used as an approximate one with a 
different type of approximate processor in order to see 
which requires greater resources. Both will be used to 
implement operators in the same one-parameter group. 
In particular, consider the operations on qudits (with or- 
thonormal basis = 1, . . . D}) specified by 



U{e) = e'%l){l\+X, 



(5.2) 



where X = J^^^^ \k){k\, and < 6* < 2tt. 

Consider the processor described by the operators Ajk 
for 1 < j,k < N, where 



_ I djkX + 6k,j+i\l){l\ J<N; 
^^■'^-i 6NkX + SkM'L){l\ j = N, 



(5.3) 



originally described in Ref. [5] . With the program state 

fc=l 



we find that for 1 < j < — 1 



1 



pi 0-1 



and for j = N 



A^(S) = -^(e^(^-i)^X + |l)(l|). 



(5.5) 



(5.6) 



What this means is that if after the action of the 
processor, the program state is measured in the basis 
{|l)p, . . . \N)p} and if the result \j)p is obtained, where 
j N, then the operation U{d) has been carried out on 
the data. However, if the result \N)p is obtained, then 
the operation U{6) has not been performed. Because 
each of these outcomes is equally likely, the probability 
of obtaining the desired result is {N — 1)/N. If instead of 
measuring the output of the program register we discard 
it, i.e. trace over it, we can use this processor as an ap- 
proximate one. The process fidelity in this case is given 
by 



F = l- 



(1 - cos{Ne)). 



(5.7) 



Another processor that will approximate this one- 
parameter group can be constructed by dividing the 
interval [0, 27r) into subintervals and approximating all 
of the operators U{9) for 6* in a particular subinterval 
by a single operator. In particular, let = ir/N, 
and approximate U{9) for 2jM < < 2{j + 1)A6' by 
Uj = U{{2j + 1)A0), where j = 0,1, . . . N - 1. We now 
define a U processor by setting, ior j,k = 0,1, ... N — 1 



Ajk — ^jkUj. 



(5.8) 



In order to approximate U{9) for 2jA9 < 9 < 2{j + 
1)A9, we choose the program state \E)p = \j)p. For this 
processor we find that 



2(i^-l). . 2(i:>-l) TT^ 

1-F< ^ ^ (l-cosA6l)^ ^ ^ 



£>2 



£)2 4iv2 ■ 



(5.9) 



By comparing the two fidelities, we see that for a fixed 
value of the program space dimension, N, the second 
processor will provide a greater accuracy. 



VI. CONCLUSION 

We have examined the approximation of a set of uni- 
tary operators by means of a programmable quantum 
circuit, i.e. a quantum processor. The programs them- 
selves are quantum states. We have shown, for a fixed 
processor, how to find the program that induces the best 
approximation of a particular unitary operator. In addi- 
tion, we have found bounds on the size of the program 
space that is necessary to approximate a set of operators 
to a given precision. 

Approximate processors can be characterized by their 
accuracy and by the resources they require. By the 
accuracy, or level of precision, we mean the quantity 
ea = 1 — min^gr niax|£5(^p) F{£,£^) [20]. Here T is the 
set of maps we want to realize, SlTip) is the set of positive 
operators on Tip with trace one (note that we are allow- 
ing mixed program states here) and £^[q\ = TvpGg^S^G'^ . 
The dimension of the program space, A^, characterizes the 
resources required. We wish to know how these two pa- 
rameters arc related. We have made some progress here 
in exploring this relation for limited sets of maps. The 
problem becomes more difficult if one considers T to be 
the set of all unitary maps and harder yet if it is the set 
of all completely positive trace-preserving maps. Once 
we have these definitions of precision and resources, we 
can consider two problems. First, given a specific degree 
of precision eq for some set of maps F, how large must 
the program space be? Second, for fixed resources, what 
is the optimal processor, i.e. for which G is the accu- 
racy the best {eg the least)? In Ref. [21] one case of this 
problem was solved by D'Ariano and Perinotti. The data 
states were qubits, and F was the set of unitary opera- 
tors acting on a single qubit. The program space was 
also a single quibit so that N = 2. They then showed 
that the optimal accuracy is given by £g = 3/4. This 
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precision can be achieved when G is a swap gate [22], i.e. 
CswAp(l^) (Xi = <E) \-ip) for all states \ip), |</)). For a 
processor with both the data and program spaces having 
the same dimension D and G given by the d-dimensional 
version of the swap gate, we find that F{U,£{), where 
£^ is the map induced on the data by the processor with 
program ^, is independent of both the program and U, 
and is equal to This implies that for this processor, 

the accuracy is given by eg = 1 — -jjj-. We suspect that 
this is the optimal value if the size of the program reg- 
ister equals to the size of the data register, i.e. N = D, 
but whether this suspicion is correct is beyond the scope 
of this paper and will be analyzed elsewhere [22]. 

There are many open issues remaining. One possibility 
is to shift our focus, and rather than ask what type of 
the processor can perform a given set of operations with 
a particular level of precision, ask instead if it is possi- 
ble to characterize the operations that a given processor 
can perform to a specified accuracy. Another issue is the 
following. So far, wc have assumed that we arc approxi- 
mating a set of unitary operators with just a single use of 
a processor. What happens if we can use the same pro- 
cessor more than once? It turns out that multiple usage 
of the processor can significantly improve the accuracy of 
the approximation. In particular, when the U processor 
(which can perform a set of unitary operators perfectly) 
is used n times, the one can perfectly perform not only 
the original set of operators, but any product of these 
operators that is of length n or less. 



It would also be useful to find specific processors, which 
arc not U processors, that can approximate a wide class 
of unitary operations. As we have seen, superpositions 
of the basis program states are not useful in optimally 
approximating a unitary operator with a U processor, 
but they very well may be useful in doing so with other 
types of processors. 

Probabilistic processors have shown themselves to be 
very flexible devices. They can perform large classes of 
operations while requiring only limited resources. Their 
drawback is that these operations are performed with a 
probability that is less than one. It remains to be seen 
how flexible deterministic processors are, but the results 
here place some constraints on what they can accom- 
plish. In this paper wc have given an example of how 
a probabilitic processor can be used as an approximate 
one. 
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